Exclusive bb pair production and irreducible background 
to the exclusive Higgs boson production 
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We calculated the cross section for the exclusive double diffractive bb production in pp — >■ p{bb)p re- 
action at LHC. Large cross sections are obtained (3 - 10 nb). This process constitutes the irreducible 
background to the exclusive Higgs production and is of particular importance in the upcoming Higgs 
boson searches. The distribution in invariant mass of the bb pair is calculated and compared with the 
corresponding contribution from the Higgs decay. The contribution from the exclusive production 
of Z° and its decay as well as the contribution from the 77 — >■ bb subprocess are also presented for 
the first time. The influence of cuts on the signal-to-background ratio is discussed. 



The identification of the Higgs boson produced in high- 
energy proton-proton colhsions is one of the main goals 
of the LHC program. The dominant mechanisms of its 
production are the gluon-gluon and WW fusions. The 
Higgs boson will be searched for in different decay chan- 
nels: bb, 77, T+r~, W'^W~ , etc., in typical inclusive 
measurements, i.e. when only the selected decay channel 
is studied, and the Higgs boson is produced in associa- 
tion with many other particles. Depending on the Higgs 
mass, the different channels seem more favorable, and 
each of the decay channels has its own difficulties for the 
experimental identification. 




An alternative solution was suggested already some 
time ago. It was proposed to identify the Higgs boson 
in exclusive process pp — > pHp in the central rapidity re- 
gion Modern estimates of the corresponding cross 
section in the fct-factorization approach with the help of 
unintegrated gluon distribution functions (UGDFs) have 
been presented in Ref . P| . It was argued that in the for- 
ward scattering limit the background in the bb channel 
should be small due to the so-called "Jz — selection 
rule" and a suppression for massless quarks in this chan- 
nel (see, for example, Ref. Q and references therein). In 
general, there is also contribution in the Jz = 2 channel 
and the 5-quarks are rather heavy in comparison with 
the typical soft hadronic scale. Therefore, a realistic es- 
timate of the background requires a real calculation for 
the genuine four-body reaction pp — > p(bb)p. The diffrac- 
tive mechanism of the exclusive bb production is shown 
in Fig. [1] For comparison, we also show the mechanism 
with intermediate Higgs boson, which is a signal in our 
analysis and a reaction with the 77 — >■ bb subprocess. 



The amplitude for the exclusive process pp — >■ p{qq)p 



FIG. 1: The diagrams for the bb pair production proceeding 
through the Higgs boson (top-left), the exclusive diffractive bb 
production (top-right) and photon-photon fusion (bottom). 
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where Ag, Ag are helicities of heavy q and q, respectively, 
ti.2 are the momentum transfers along each proton line, 
qi,t,q2,t,xi^2 and (70, t, x'l ^ x'2 xi^2 are the transverse 
momenta and the longitudinal momentum fractions for 
active and screening gluons, respectively. Above fg^i/2 
are the off-diagonal UGDFs of both nucleons. In the cal- 
culations presented here we take /x^ = M^j^, where M^i 
is the invariant mass of the bb system. This is consis- 
tent with using //^ — Mfj — M^^ for exclusive Higgs 
boson production, which is preferred from the theoret- 
ical point of view [3]. The hard g*g* — > bb subpro- 
cess amplitude Va^a^ consists of the Higgs decay signal 
9*9* H ^ bb and background contributions, where 
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the major ones come from the direct g*g* — bb process, 
from the photon-photon fusion and from the exclusive 
production of and its subsequent decay Z° — > bb. 
Note, that due to integration over qg ^ in the diffractivc 
amphtude ([1}, only symmetric part of the hard subpro- 
cess amplitude VA^A,-(qo.i) = Vx^x-{-c[o^t) contributes to 
the diffractive cross section. 

The off-diagonal UGDFs are written as 

/■off/ / 2 22\D-f/ 2 2\ /r,\ 

fg [x ,xi^2, qi^2t^<lQ,t^fJ'F) - Rg Jgixi.2,qi,2t^^J-F)i (2) 

where Rg ~ 1.2 which accounts for the single logQ^ 
skewed effect In the considered kinematics the di- 
agonal unintegrated densities can be written in terms of 
the conventional (integrated) densities xg(x,qf) as 
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[xgix,q^)^/U^)] , (3) 



where Tg is the conventional Sudakov survival factor 
which suppresses real emissions from the active gluon 
during the evolution, so the rapidity gaps survive. The 
gluon qt 's typical for the central exclusive Higgs produc- 
tion at LHC are of the order of few GeV 0] ■ 

Additionally, following to Ref. we use the factor- 
ization scale 11 F — J^bb t given by the transverse mass 
of the bb pair My^ t compared to the KMR conven- 
tion Q /x™^ = M^i j/2. We will discuss uncertainties 
due to sensitivity to the factorization scale /i p choice and 
conventional gluon densities xg{x,q^) at small x and qt 
below when presenting numerical results. 

In the framework of the /cj-factorization approach [lol | 
the hard subprocess g*g* — >■ qq gauge invariant ampli- 
tude reads 



V^X{qi.q2) ^ - ^ 
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where Ep cms — \/s/2 is the c.m.s. proton energy, are 
the color group generators in the fundamental represen- 
tation, u{ki) and v(k2) are on-shell quark and antiquark 
spinors, respectively, &, b are the effective vertices arising 
from the Feynman rules in quasi-multi-Regge kinematics 
(QMRK) approach iUustrated in Fig. [2j 
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where r'"^'^((ji, (J2) is the effective three-gluon vertex. 
These effective vertices were initially proposed for mass- 
less quarks in Ref. [ll| and then extended for massive 



case in Ref. 13, 12 1. The effective ggg- vertices are can- 
celed out when projecting the qq production amplitude 
Eq. (jH) onto the color singlet state, so only the first two 
diagrams in Fig. [2] contribute to the final result for the 
production amplitude. Since we will adopt the definition 
of gluon polarization vectors proportional to transverse 
momenta qi/2±, i-e. £1,2 ~ '?i/2±/3^i,2 (see below), then 
we must take into account the longitudinal momenta in 
the numerators of effective vertices (see Eq. ([5])). 

The SU(3) Clebsch-Gordan coefficient (3i,3fc|l) = 
5"^^ / y/Wc in Eq. (U) projects out the color quantum num- 
bers of the qq pair onto the color singlet state. Fac- 
tor 1 / yJNc provides the averaging of the matrix element 
squared over intermediate color states of quarks. 

Therefore, we have the following amplitude 
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Taking into account momentum conservation and us- 
ing the gauge invariance property, we get the following 
projection to the light cone vectors (so called "Gribov's 
trick") 
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Last expression shows that an important consequence of 
the gauge invariance is the vanishing of the matrix ele- 
ment of the effective ggqq-veTtex between on-mass-shell 
quark and antiquark states in the limit of small qi± and 
q2± [il[l2| 



for 



qii_ or q2A 



0, 
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The amplitude ([6]), projected out onto a particular 
quarkonium state, was successfully applied for the de- 
scription of the recent CDF data on the exclusive produc- 
tion of charmonia in Ref. [31 ■ Here we apply the same 
formalism for separate b and b jets production. This is 
exactly the same formalism as used recently for the ex- 
clusive open charm production in Ref. [6] . 

In the experimentally important case of the forward 
proton scattering p'^ 2t ~^ we have qi,t — — q2,t — 
qo,t = qt, and the transverse mass of the bb pair in terms 
of 6-quark rapidities ^ reads 



Mb-b,t = ^4 = ,(1 + cosh(yb - Vb)) 



(8) 



If one looks at centrally produced 6-jets only, i.e. j/^ — >■ 
0, then according to Eq. ([8|) the only way to produce the 
large invariant mass M^j, ~ Mh is to consider high-fcj jets 
limit kt ^ mb,qt- The amplitude vanishes for jets 
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FIG. 2: Effective vertex in the QMRK approach [ffj]. Last diagram with effective 3-gluon vertex drops out in projection to the 
color singlet final state. 



with very small rapidities, i.e. when yb ^ Vi ^ 0. Am- 
plitude V++ behaves as ~ cos (j)/k^ and thus extremely 
suppressed in the high-Zcj limit. This is in agreement with 
the statement that Higgs CEP background is suppressed 
in very forward and quark massless limits for centrally 
produced hh jets (with small rapidities) , and agrees with 
the Jz = selection rule [15 1. 



However, the particular high-fct limit is only a part of 
the whole story. In our previous analysis of the exclusive 
open charm production in Ref. @ it was shown that the 
dominant contribution to the cc dijet cross section comes 
from relatively small quark transverse momenta fcj ~ 1 
GeV. The same should hold for exclusive hh pair produc- 
tion relevant for the Higgs background. Indeed, resolv- 
ing relation ([5]) with respect to typical rapidity difference 
\yb — yb\ = ^y neglecting quark transverse momenta kt 
and keeping only the 6-mass contributions rrib ~ 4.5 GeV 
at fixed M^^ = 120 GeV, we get ~ 6.6. So, the ir- 
reducible background for Higgs CEP can be dominated 
by 6-quarks with comparably small transverse momenta 
kt <^ rrib, but with rather large rapidities yb ^ —yi — 3.3. 

In the last kinematical situation the amplitudes Va^a^ 
are not suppressed by a large denominators, and signifi- 
cant contributions can be obtained. For simplicity, keep- 
ing only the quark mass ^ kt and large rapidity 
Vb ~yb ^ contributions, we get for helicity ampli- 
tudes 
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where we have introduced the shorthand complex nota- 
tions in the forward limit 

K^ky+ik'^, p = qy + iq^^ 

for the quark and gluon transverse momenta, respec- 
tively. 

In this \ow-kt regime helicity amplitude is domi- 
nated over V++ as opposite to the high-fcj case. From 
Eq. ^ we see that the amplitude is proportional to 



qt = — IpP, which typically can be of the order of few 
GeV at LHC energy. This means that numerically low- 
kt contribution (of course, at not extremely large y^i) in 
the case of 5-jets can lead to a dominant contribution to 
the exclusive background for Higgs CEP. In this asymp- 
totics, the quark mass rUb plays an important role since it 
comes into the denominator in Eq. (jO)). Precise evaluation 
of the corresponding signal, however, demands employ- 
ing the formulae for the hard amplitudes in the general 
form in Eq. More detailed analytical and numerical 
investigation of contributions from different parts of the 
phase space will be presented elsewhere (l6| . 

In parallel to the total cross section, we calculate the 
differential cross sections for exclusive Higgs boson pro- 
duction. Compared to the standard KMR approach here 
we calculate the amplitude with the hard subprocess 
9*9* ~^ H taking into account off-shellness of the active 
gluons, i.e. fully consistent with the exclusive production 
of the hh pairs, where the gluon transverse momenta play 
crucial role. The details of the off-shell matrix element 
can be found in Ref. [l^ • In contrast to the exclusive Xc 
production [3] , due to a large factorization scale ~ Mh 
the off-shell effects for g*g* ^ H give only a few percents 
to the final result. 

The same unintegrated gluon distributions based on 
the collinear distributions are used for the Higgs and 
continuum hb production. This is absolutely necessary 
for proper estimate of the signal-to-background ratio, 
the main purpose of the present Letter. In the case 
of exclusive Higgs production we calculate the four- 
dimensional distribution in the standard kinematical 
variables: y,ti,t2 and 4>- Assuming for this presentation 
the full coverage for outgoing protons^ we construct the 
two-dimensional distributions da/dycPpt in Higgs rapid- 
ity and transverse momentum. The distribution is used 
then in a simple Monte Carlo code which includes the 
Higgs boson decay into the 66 channel. It is checked sub- 



^ The exact comparison with experimental set up requires inclu- 
sion of extra cuts on fractional momentum loss of protons and is 
shghtly different for ATLAS and CMS. 
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FIG. 3: The bb invariant mass distribution for ^/s = 14 TeV 
and for -2.5 < y < 2.5 corresponding to the ATLAS/CMS 
detectors. The absorption effects were taken into account by 
multiplying by the gap survival factor So = 0.03. 
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FIG. 4: Illustration of the gluon densities as functions of the 
longitudinal momentum fraction x at the characteristic fac- 
torisation scale = 2 GeV^ (typical for the central exclu- 
sive production processes) given by the global parameteri- 
zations CTEQ6L1 [H], GRV94HO [H, MSTW2008LO and 
NLO M. 



sequently whether the b and b enter into the pseudorapid- 
ity region spanned by the central detector. Including the 
simple cuts we construct several differential distributions 
in different kinematical variables. 

In general, employing the diagonal UGDF in the form 
([3]) one encounters a problem of poorly known gluon 
PDFs at rather low xi,2 and especially small gluon virtu- 
alities q'^ . For an illustration of the corresponding uncer- 
tainties, in Fig. |4]we show several parameterizations for 
the gluon PDFs widely used in the literature as functions 
of momentum fraction x at the evolution scale ~ fixed 
at characteristic value 2 GeV^ typical for the exclusive 
production of Higgs boson (for quarkonia production it 
is even smaller leading to huge uncertainties as discussed 
e.g. in Ref. We see that at x < 10"^ the PDF 



uncertainties may strongly affect predictions for not suf- 
ficiently large gluon transverse momenta. In this sense, 
precise data on the diffractive and central exclusive pro- 
duction can be used to constrain the PDF parameteriza- 
tions 
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Testing other models of UGDFs different from Eq. ([3]) 
may be important for estimation of an overall theoretical 
uncertainty of our predictions and their stability, and it 
is planned for our future study. 

In the following we shall present the main results. A 
more detailed analysis with the presentation of several 
differential distributions will be given elsewhere [H. In 
Fig. [3] we show the most essential distribution in the in- 
variant mass of the centrally produced bb pair, which is 
also being the missing mass of the two outgoing protons. 
In this calculation we have taken into account typical de- 
tector limitations in rapidity —2.5 < yb,yb < 2.5. We 
show results with different coUinear gluon distributions 
from the literature: GRV CTEQ fl^, GJR 21] and 
MSTW [13. The results obtained with radiatively gen- 



erated gluon distributions (GRV, GJR) allow to use low 
values of Qt — qo,t, <}i,t, <l2,t whereas for other gluon dis- 
tributions an upper cut on Qt is necessary. The inte- 
grated double-diffractive bb contribution calculated here 
seems bigger than the contribution of the exclusive pho- 
toproduction of bb estimated in ;22] and details require 
systematic studies in the future. The lowest curve in 
Fig. |3] represents the 77 contribution (the bottom dia- 
gram in Fig. [ij. While the integrated over phase space 
77 contribution is rather small, is significant compared 
to the double-diffractive component at large M^f, > 100 
GeV. This can be understood by damping of the dou- 
ble diffractive component at large M^f, by the Sudakov 
form factor [1, 0| • In addition, in contrast to the double- 
diffractive component the absorption for the 77 compo- 
nent is very small and in practice can be neglected. 

In the top panel of Fig. [5] we show the double diffractive 
contribution for a selected (GTEQ6 [3]) collinear gluon 
distribution and the contribution from the decay of the 
Higgs boson including natural decay width calculated as 
in Ref. [H], see the sharp peak at M^j = 120 GeV (as- 
sumed arbitrarily for illustration), which is not excluded 
at present by the Higgs searches at LEP and Teva- 
tron [2^. The phase space integrated cross section for 
the Higgs production, including absorption effects with 
Sg = 0.03 is slightly less than 1 fb which is similar to that 
predicted by the KMR group. This value is similar as in 
many KMR evaluations [4]. The result shown in Fig l5] in- 
cludes also the branching fraction for BR(77 — >■ bb) ~ 0.8 
and the rapidity restrictions. The second much broader 
Breit-Wigner type peak corresponds to the exclusive pro- 
duction of the Z° boson with the cross section calculated 
as in Ref. [2^. The exclusive cross section for ^/s = 14 
TeV is 16.61 fb including absorption (28.71 fb without ab- 
sorption effects). The branching fraction BR(Z° — ^ bb) sa 
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FIG. 5: The bb invariant mass distribution for ^/s = 14 TeV 
and for b and b jets in the rapidity interval —2.5 < y < 2.5 
corresponding to the ATLAS detector. The absorption effects 
for the Higgs boson and the background were taken into ac- 
count by multiplying by the gap survival factor So = 0.03. 
The top panel shows purely theoretical predictions, while the 
bottom panel includes experimental effects due to experimen- 
tal uncertainty in invariant mass measurement. 



0.15 has been included in addition. In contrast to the 
Higgs case the absorption effects for the production 
are much smaller [26| . The sharp peak corresponding to 
the Higgs boson clearly sticks above the background. In 
the above calculations we have assumed an ideal no-error 
measurement. 

In reality the situation is, however, much worse as both 
protons and in particular b and b jets are measured with 
a certain precision which automatically leads to a smear- 
ing in MjjI . While such a smearing is negligible for the 
background, it leads to a significant modification of the 
Breit-Wigner peaks, especially of the sharp one for the 
Higgs boson. In the present Letter the experimental ef- 
fects are included in the simplest way by a convolution of 
the theoretical distributions with the Gaussian smearing 
function 



G{M) 



1 



(10) 



with (7 = 2 GeV, which reahstically represents the ex- 
perimental situation 27|, |28| and is determined mainly 
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FIG. 6: The bb invariant mass distribution for ^/s = 14 TeV 
for a limited range of 6 and b rapidities: — 1 < y < 1 (left) 
and for pi,2t < 0.4 GeV (right). 



by the precision of measuring forward protons. In the 
bottom panel we show the invariant mass distribution 
when the invariant mass smearing is included. Now the 
bump corresponding to the Higgs boson is below the bb 
background. With the experimental resolution assumed 
above the identification of the Standard Model Higgs will 
be rather difficult. The situation for some scenarios be- 
yond the Standard Model may be better [l6| . 

The question now is whether the situation can be im- 
proved by imposing further cuts. In Fig. [5] (left panel) we 
show the result for a more limited range of b and b rapid- 
ity, i.e. not making use of the whole coverage of the main 
LHC detectors. Here we omit the Z° contribution and 
concentrate solely on the Higgs signal. Now the signal- 
to-background ratio is somewhat improved. This would 
be obviously at the expense of a deteriorated statistics. 
Similar improvements of the signal-to-background ratio 
can be obtained by limiting transverse momenta of out- 
going protons (right panel). 



117.6< M^^< 122.4 GeV 



PP— >pHp, H— >bb M„ = 120G8V 




FIG. 7: Two-dimensional distributions over b and b rapidities 
for QCD background (left panel) , and Higgs CEP signal (right 
panel) with the optimal two-dimensional cut in the {yt, yi) 
space is marked by the thick contours. 

In order to preserve statistics and to remove most of 
the bb background we have to impose more specific two- 
dimentional cuts. An example is shown in Fig. [T] Indeed, 
considering only the band between two thick solid lines 
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would remove most of the bb background, concentrated 
mainly in regions with relatively large difference between 
quark and antiquark rapidities \yb — Vil ^ 1 (see, the left 
panel in Fig. [7]). At the same time, such a cut allows 
to keep most of the Higgs signal, which is concentrated 
in the central rapidity region j/f, ^ « 0, as oppose to the 
background (see, the right panel in Fig. [7]). 

In the present analysis we have not been interested 
in the precise estimation of the cross section but rather 
in understanding the signal-to-background ratio which is 
of the major importance for the upcoming Higgs boson 
searches at the LHC. Consequently, we have presented 
results with only one UGDF. This ratio is practically the 
same for other UGDFs, which will be shown explicitly in 
[lij . The absorption effects have been included here in a 
simple multiplicative form. They are expected to be the 
same both for the signal and the background, and thus 
are not affecting the ratio under consideration. The same 
gap survival factor has been used in both cases. 

The overall theoretical uncertainty of our predictions 
for the absolute value of the bb background contribution 
is estimated to be the same as for the Higgs CEP, and 
given by a factor of 3 more/less Theoretical uncer- 
tainty in the signal-to-background ratio under consider- 
ation is typically much smaller, as the main part of the 
uncertainty coming from the normalisation of UGDFs 
{Rg factor) and absorbtive effects (gap survival Sq fac- 
tor) is common for both contributions and thus canceled 
out in the ratio. 

In our analysis we have been concentrated on the 
irreducible background only. Other contributions, al- 
though, in principle, reducible, can in practice be also 
rather troublesome [29|. These include dijet misidentifi- 
cation (mainly due to the qq — » gg subprocess), inclusive 



double-pomeron processes 
lated to large luminosity 



and multi-event effects re- 
. Further analyses, espe- 
cially for the Standard Model Higgs boson production, 
seem to be necessary to understand whether the Higgs 
boson can be identified in the exclusive production, per- 
haps not only in the bb decay channel. The present par- 
ton level analysis should be supplemented in the future 
by additional analysis of bb jets by including a model of 
hadronization. Then standard jet algorithms could be 
imposed and the quality of the b and b kinematical re- 
construction could be studied in detail. 

We are indebted to Valery Khoze, Oleg Teryaev for 
a discussion of central diffraction, to Andy Pilkington 
and Christophe Royon for a discussion of experimental 
limitations of the exclusive processes discussed in this 
letter and Ania Cisek for providing us numerical results 
of the exclusive production from Ref. [26| . This work 
was partially supported by the Polish grant of MNiSW 
N N202 249235 and by the Carl Trygger Foundation. 
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